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Integer filling factors of a hetero-structure 2D electron gas (2DEG) in 
an external magnetic field are special ones because a huge degeneracy of 
the ground state is gone here. It justifies the Hartree-Fock approximation 
with the accuracy limited only by normally a small parameter: V^''^* /huo, 
where V^^^ is the energy of the Coulomb interaction and uo is the frequency 
of the cyclotron resonance. Such an approach predicts the ground state 
of a single layer 2DEG at = 1 to be a ferromagnet with the degenerate 
total spin orientation. The elementary excitations of 2DEG are electron-hole 
pairs or excitons and in the close binding limit of vanishing momentum they 
transform into the elementary excitations of a ferromagnet - spin-waves. The 
latter are gapless [|l| and do not interact with each other p| if Zeeman energy 
is neglected - the two consequences of the Goldstone theorem. In the limit of 
large momentum the electron and the hole of an exciton are well separated 
and they become the elementary charged excitations. In the Ref.0 it was 
shown that a special topological spin texture in a 2D ferromagnet called a 
skyrmion |Q carries a unit charge while costing only half the exciton's energy 
to appear. 

The case of a bilayer 2DEG turned out to be a more rich one. The Hartree- 
Fock approximation does not apply here except for two limiting cases. The 
first one is the case of well separated layers which is a common setup in the 
experiment |^ and where, theoretically, one starts from the two single layer 
ferromagnets and makes the perturbation expansion in powers of interlayer 
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interactions [0. And the second one is the symmetric case defined in such 
a way that one can freely rotates an electron spinor in both layer and spin 
spaces. The latter requires to approximate the Coulomb interaction by its 
symmetric part and to neglect all symmetry breaking fields like Zeeman en- 
ergy. The first attempts in this direction dealt with the case of filling factor 
u = 1 and relied heavily on the assumption of a saturated spin polarization 
of electrons ||, M. This symmetric approximation turned out to be useful to 



determine the exciton energy in bilayer M. Recent works Refs.jlO, 11| spe- 



cialize to the bilayer heterostructure case u = 2, employ the Hartree- Fock 
approximation and predict a phase diagram that features three phases: the 
ferromagnetic, the canted antiferromagnetic and a special spin-singlet phase. 



In this paper we reproduce the phase diagram of the Refs. |TD|, |TT| isolating the 
symmetric and the symmetry breaking parts of the Hamiltonian in a consis- 
tent way. Our approach reveals the Hartree-Fock phase diagram to be indeed 
exact in the limit V"-™'^ /V^y"^ 0, where y?'™- is the 5'f/(4)-symmetric part 
of the bilayer Hamiltonian whereas ^Z"^"** is anisotropy interactions that re- 
duce the bilayer Hamiltonian symmetry to SU{2) ® SU{2). We prove the 
stabihty of all phases with respect to long-range spatial perturbations. 

Our new results concern the energy of topological excitations in bilay- 
ers. We find that low-energy excitations over the bilayer ground state is 
governed by the f/(4)/?7(z/) ® ?7(4 — u) coset non-linear Sigma Model. We 
identify charged excitations in the bilayer with skyrmions or topological exci- 
tations of this non-linear Sigma Model i, I, m. We calculate the skyrmion 
energy gap to vary dramatically over the bilayer parameter space and we 
find a sharp dip of this gap in the canted antiferromagnetic phase of the 
bilayer. Our work was motivated by recent measurements of the diagonal 
conductivity activation energy 0. In this paper we suggest to identify the 
experimental activation energy with the energy gap of skyrmion, and we have 
got a qualitative agreements with the results of Ref.[B] along this line. 



Hamiltonian of 2DEG bilayer 

The electronic Hamiltonain of a 2DEG in a confining potential V{p) and 
in an external magnetic field H consists of a one-particle part as well as a 
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Figure 1: Schematic view of a confining potential V{^) in a typical Bilayer 
setup 



Coulomb interaction part: 

(1) 

where a, (3 — ± are spin indices and thereafter a sum over repeated indices 
is imphed. We use such units that h — 1, e — c and H — B — 1. The latter 
implies that all distances can be expressed in terms of the so-called magnetic 
length: Ih — cfi/ eH = 1. We split three coordinates p into a perpendicular 
to the layer coordinate ^ and two in-plane coordinates r = {x, y) = [z, z). We 
assume that the confining potential is uniform over the plane: V{p) = V{^), 
and represents a double well structure in the transverse direction as shown on 
the Fig.l, with the two wells being separated by the distance d. We use only 
two eigen functions: the lowest energy symmetric XsiO ^^^^ antisymmetric 
XaIOj fi'oni a set of one-electron eigen functions in the confining potential 
V{^) and we expand an electron second-quantized operator in terms of these 
two eigen functions: 
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(2) 



T,n,p 



where c^^^ and c^rp are electron creation and anhilation operators, (j)n,p{zz) 
is an electron wave function number p in the Landau gauge in the n's Landau 
Level, the index r = 1, 2 being the layer index and the layer wave functions 
read: 

xiAO = ^ • (3) 

We restrict our model to the case of a sufficiently strong magnetic field, such 
that the cyclotron energy huo dominates over the Coulomb, Zeeman and the 
level splitting: Ea — Es, energies. Thus, we specialize to the lowest Landau 
Level and retain only the term n = in (H). 

Plugging the wave function into the Hamiltonian (0) we find a 2DEG 
Hamiltonian as: 
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P1--P4 



^aripi ^I3t2P2 ^P^SPa ^^Q-taPa i 



where we have defined a hopping constant: 



(4) 



^ = \j j d'rd^ ^li^XrAiXr^ViOXr^mpi^. (5) 

as well as an external electrostatic potential created by an asymmetric gate 
charge: 

whereas the Coulomb interaction matrix reads: 



We use notations: r^, and r^, for the Pauli matrices in the layer space 
whereas we use notations: o"^, and cx^, for the Pauli matrices in the spin 
space. The hopping constant can be related to the splitting of the symmetric 
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and the antisymmetric levels: t = Ea — Es- The electrostatic potential /i^, 
which can be viewed as a difference between the chemical potentials in the two 
layers, breaks down the symmetry between the two wells of V{S,) potential. 
This term appears naturally when a single gate is fabricated to control the 
electron density in the bilayer. In the limit d ^ 0, fi^ vanishes too, whereas in 
the limit of large layer separation: d — >■ oo, /i^ — > cxo and electrons reside only 
on the layer adjacent to the gate. We assume that the energy of a capacitor 
formed by the two layers is much lower than the characteristic Coulomb 
energy e'^/nl^, per area, where k is the dielectric constant. We note the 
invariance of the Coulomb energy: (|^ under the following transformations: 
Ti r4, T2 ^ T3 as well as {tiT4^) ^ {T2T3). To fully exploit these symmetries 
we cast the Eq.|^ into a more suitable representation: 

where r° is the unit matrix, V^'^ is a 3 x 3 symmetric interaction matrix with 
indices /i, v running over a set (0, z, x). If there is a symmetry of the Coulomb 
interaction under an exchange of layers: {^E,') ^ (— ^ — and 1 2 then 
it restricts further values of the interaction matrix: V^^ = and V^^ = 0. 
But in the presence of a gate asymmetry we shall keep the matrix element: 
V^^. Therefore, the Coulomb interaction matrix for symmetric bilayer 2DEG 
depends on four parameters: V^°° > 0, V'^^, V^^ > 0, V^^ > 0. We note also, 
that V^'^ ~ X) ~ X^) whereas V^^ ~ d'^/\z\^ as \z\ 00. In the following 
we shall neglect V^^ matrix element. 

Next, we split the total bilayer Hamiltonian (^) into two parts: the first 
one contains a dominant Coulomb energy term: 

H''"" = ^ct,c^r, + IJ ^V'\q)N{q)N{-q), (9) 

where V^'^{q) is the Fourier transform of V^^ {r) multiplied by a factor 
exp(— and the electron density operator reads: 

N{q) = J2 ^trpCarp-qyeW-iqxip-'^) (lO) 

p ^ 

This part of the Hamiltonian is invariant under uniform rotations from the 
SU{A) Lee group in the combined spin and layer space. Every of its eigen 
energy is hugely degenerate. Given any eigen state |\E')o a set of related 
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eigen states can be generated by applying rotations: |\E') = U\'^)o, where 
U e SU{4). For Landau level filling factor u = 1, u = 2 and u = 3 we 
assume that the the bilayer ground state is uniform over ]9-orbitals: 

1/ 

^ = niKr.p \empty), (11) 

i=l V 

and we prove in the next section that this state is stable with respect to long- 
range spatial perturbations. One can easily check by inspection that any 
such wave-function (pJ]) represents an eigen function of the if**^™ (^. The 
remaining few terms in the Hamiltonian @) are treated like perturbations: 

+ \j ^.v'^^m^'m'-i-q), ^^^^ 

where (see e.g. 0) 

T'iq) = E c+,,X,^c„.,,_,,exp-zg.(p-|) (13) 

with {fiu) 7^ (00). The Hamiltonian (|T2|) breaks down the SU{4) symmetry 
but it is still invariant under separate rotations in the spin and layer space: 
SU{2) ® SU{2). We shall call this part of the Hamiltonian the anisotropy 
Hamiltonian. It lifts the degeneracy of eigen states of the SU (4)-symmetric 
Hamiltonian (|^). An important point to note here is that a splitting of energy 
levels is determined by matrix elements of weak anisotropy Hamiltonian (0) 
truncated to a linear space of the symmetric Hamiltonian (^ level degener- 
acy. There are no Fermi-liquid type renormalizations of the constants of the 
anisotropy Hamiltonian (0) due to the S'f/(4)-symmetric Hamiltonian (|^). 
In other word the mean- field Hartree-Fock approach is perfect for the = 1, 
u = 2 and u = 3 cases. 

Our guiding analogy in treating the total bilayer Hamiltonian (^, |l^) 
lies in the theory of magnetism. We will see below that there exists a local 
order parameter: Q, very much like magnetization. And we aim to express 
the total bilayer Hamiltonian (§,0) in terms of this order parameter Q. 
The exchange-like Hamiltonian (^ has to be expanded in powers of spatial 
variations of order parameter Q{r) with the second power of gradients being 
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the important contribution, whereas only locally homogeneous Q has to be 
retain in the anisotropy Hamiltonian ([121) . In the next section we carry out 
the first step whereas in the next-to-next section we transform the anisotropy 
energy. 



SU(4) Symmetric Case 



In this section we specialize to the SU (4)-symmetric part of the bilayer 2DEG 
Hamiltonian (^ which is invariant under the global rotations of a four com- 
ponent electron spinor by the 4x4 matrix U from the 5*^7(4) Lee group. 
We find it useful the ground state of operators c to be the reference state. 
Any non- homogeneous state is then generated by rotation: U{t,r). And the 
action of bilayer is some functional of it: 

S[U{t,r)] = -i trlog J Vc+{t)Vc{t) exp (^i J Cd?j (14) 

where the symmetric Lagrangian of bilayer 2DEG reads: 
d 



(fq 



C/3r2 d r + 



(15) 



V''{q}N{q}N{-q}, 



where Ao = (0, A^) is the vector potential of the external magnetic field and 



d 

n' = zu+—u, Q 

ot 



-I u+vu. 



(16) 



In the limit of slow spatial variations of rotation U{t, r) this functional has 



an expansion in powers of Vt. In this context the functional (|l^) is called 
an effective low-energy Goldstone Action, and we are going to find it in this 



section. All calculations follows step in step those done in Ref. fl^ for the 
case of a single layer, and here we emphasized only the points of difference. 

Calculations of the Goldstone Action can be carried through for three 
filling factors of the bilayer 2DEG: z/ = 1, z/ = 2 and z/ = 3, at once. We 
start with choosing the reference state: i) = 1, one electron with spin up 
fills every orbital p of the lowest Landau Level of the first layer: ii) 



7 



n = 2, two electrons - one with spin up on the first layer and the other with 
spin down on the second layer - fill every orbital p the lowest Landau Level: 
|+)i|~)2; iii) three electrons - one with spin down on the first layer and the 
two other with spin up and down on the second layer - fill every orbital p 
the lowest Landau Level: |— )i|+)2|— )2- The case iii) reduces to the case i) 
if one makes the electron-hole transformation. 

The one-electron Green functions defined for the reference state of the 



Hamiltonian ([15D in the homogeneous limit: = 0, reads: 



C«/3,TlT2 



-'a/3,Tir2 



— 

P,T\T2 al3,TlT2 



e + Eo 



iO 



iO' 



where /z is the chemical potential, 



and in the case n 



o/3,Tir2 



1,3: 



^ 1 

10 

10 

V 1 y 



''afS,TiT2 



whereas in the case n = 2: 



± 



/I 




V 
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-1 / 



-'a/3,TiT2 
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The effective action can be written as: 

S* = 5*0 + 5*2, 

where 



So = iti log 



G_ 



(17) 



(18) 



(19) 



(20) 



(21) 
(22) 
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Figure 2: Second order Hartree-Fock diagrams 



and 5*2 is conveniently represented in terms of diagrams in the Fig.2. As it is 
explained in the Ref. [T^ the first order perturbation correction to the Green 
function: 5G = G — Gq, contains an electron propagation on the first excited 
Landau Level. Thus, we also need the bare Green function of an excited 
electron on this level: 



2 



/3,T1T2 



Ol(S,TlT2 



1 



■'al3,TiT2 



e — 1/m + El — /i + iO 



where 



Eq = 2Ei 



e — 1/m — fi + iO' 



71 e 



2 kI 



(23) 



(24) 



H 



As it is explained in the Ref. [Q the gradient vector field enters the following 
terms in the Hamiltonian 



^ 2m 



i (^m^m) 



<^f^,TlT2 



Cf3T2 d'^r, 



(25) 
(26) 
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where the operators Il± shift an electron between the adjacent Landau Levels: 



n_ 



^2ra0„_ip(f) 



n 



(27) 



though only the n = and n = 1 Landau Level states are relevant for our 
problem. The two new gradient vector fields in the Hamiltonian (pSf ) are 
defined as follows: 

n± = -tU+{dyTtd.)U, (28) 

These two components are real: f2_ = An expansion of the 2DEG action 
up to the second power of the Hamiltonian ( p5D reads: 



SSo = tti (HiGo) + -tr (HiGoHiGo) + tti {H2G0) 



(29) 



Combining this expansion with diagrams of the Fig. 2 we find the low-energy 
Goldstone Hamiltonian as follows: 



Hg = yJ |^tr((S°-S^)l]_(r-)(S° + S^)f]H 



(30) 



The insertion matrices in ( pOD are non-negative diagonal ones and they rep- 
resent the occupation number for the electron states: 



N = - 
2 



and 



N 



9 V 






















(31) 



(32) 



where blocks are 1x1 and 3 x 3 in the case z/ = 1, 3 and 2 x 2 in the case 
u = 2. It follows immediately that Hq > 0. 

The matrices N and 1 — can be viewed as projector operators that allow 
only those rotations into the Hamiltonian (^) that do change the ground 
state. It is useful to understand what particular sub-set of SU (4) Lee group 
these physical rotations form. The vector field can be expanded in the 
basis of fifteen generators of 5'f/(4) Lee group: {S'}, with I = 1..15. And 
we subdivide them into two complementary sets: the first one includes those 
generators that do commute with occupation number matrices ( |3ll , ^2|) , and 
we shall called it an [even) set, whereas the second one includes the remaining 
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generators, and we shall call it the (odd) set. Generators of the (even) set 
constitute an algebra itself. This algebra has a normal Abelian subalgebra 
formed by a single generator: — trS^/4. A Lee group built around the 
(even) set of generators is called a stabilizer sub-group 5" of SU (4) Lee group. 
The (odd) set always contains an even number of generators. Specifically, 
eight in the case of z/ = 2 and six in the case u = 1,3. 

The Hamiltonian (^) must be invariant under the time reversal sym- 
metry. The time reversal operator can be chosen as a complex conjugate 
operator: U — » f/*. It follows that Q —Q'^. Now it is evident that f2_ and 

get interchanged under the time reversal in (^), and this changes the 
Hamiltonian. But we remember that the time reversal is always accompanied 
by inverting the magnetic field and, thus, we can restore the time reversal 
symmetry by multiplying the antisymmetric part of the Hamiltonian Hq by 
the sign of the magnetic field: 

i^G = Y / ^ [tr - N)n^Nn^) + I sgn(5^)e^,tr (fi^fi.AT)] (33) 

Now if we define a useful gradient vector field: 

= ^tr {u^m^d,U{r)) , (34) 

then it is straightforward to rewrite the Eq.(p3D: 

^G = Yj^b' ((S° - N)n^Nn^) + sgn(5^) curl fi^] , (35) 

The first term in Hq (|35| ) is the gradient energy whereas the second term is 
proportional to the topological index of an excited state: 

Q = [ curlfi^ — = Z, (36) 

where Z is the set of integer numbers. The case Q = ±1 corresponds to the 
simplest spin skyrmion in the first layer being rotated by a SU{4) matrix 
to become a general bilayer skjnrmion. The energy constant in Hq (^) 



coincide identically with that of the one-layer case |]T2|, which means that 
the bilayer skyrmion energy is the same as found for one layer. But there 
is an important difference between a bilayer skyrmion and a spin skyrmion 
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in a single ferromagnetic layer. Namely, as was shown in the Ref.[|T2| any 
skyrmion carries a charge density in the core: 

n(r) = = 7j, (37 

^ ' 2n 7r(i?2 + 7=^)2' ^ ^ 

where R is the radius of the skyrmion's core, contrary to that in the bilayer 
a charge of the skyrmion's core can be delocalized over the two layers with a 
much longer tails: 

The total charge of two layers does converge according to (^) at large dis- 
tances from the skyrmion's core. 

The Goldstone Hamiltonian (^) can be cast in a special order parameter 
representation. To do this we define a non-homogeneous order parameter 
matrix Q as follows: 

Q(f) = U{r)NU+{r). (39) 
This electronic order parameter has an important property: 

{A)=tTiAQ), (40) 



where A is any operator. Inspecting the particular difinition of N ( plD it 
becomes evident that rotations from the denominator sub-group S leaves the 
order parameter intact. Thus, rotations in (pQ]) can be restricted to a coset 
or, in other word, a physical space of the bilayer 2DEG: 

m . (41) 

U{u) ® f/(4 - z/) 

Now it a straightforward calculation to rewrite Hq (^Sj) in terms of the order 
parameter matrix: 

H=^JtT (Vgvg) ^ + sgn(5^)^ J (Qd.Qd^Q) ^. (42) 

In this representation the topological index appears as an index of a map of 
the order parameter coset space into a 2D plane. The index selection rule 
( |36D is a consequence of a well known homotopy group identity: 
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In the end we have to include the Coulomb energy of charge distribution 
inside the skyrmion core: 



irrr ^ f /" j2 -j2 ^, curl fi^ (r) 6^ curlfi'^(r*) 

^Ha = - J J Mr' (44) 

Anisotropic Part of Coulomb Energy. Phase 
Diagram. 

In this section we cast the anisotropic part of the bilayer Hamiltonian (0) 
in terms of the order parameter matrix Q. It can be conveniently done by 
the following Harteree-Fock average of c-operator product in (|12|): 



T1T4 T2 

where r'^ acts on four-spinor as r'^ ® cr^- Next, we define the following 
Coulomb anisotropy constants: 

where the last approximation holds for (ab) ^ (00) in the limit d -C Ih- 
And, finally, we rewrite the anisotropy Hamiltonian ([T2|) in terms of order 
parameter matrix Q: 

-\g\piBHii {Qa') + ^^^^ [tr {Qr') tr (Qr^) - tr {Qt'Qt')] , 

(47) 

where M is the number of degeneracy of the Landau Level. The Eqs.(^ 
defines the effective long-range Hamiltonian of a bilayer at integer filling 
factors. At non-zero temperatures thermal fluctuations of the order param- 
eter soften the anisotropy constants in the Hamiltonian (|47|) . The relevant 
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calculation can be found in eg. |T3[ and the result reads: 




8T/_Ei 




V 



V 



R 



R 




V — 



V — 



8T/_Bi 



(48) 



{\9\iJ'bH)r 



( 



R* 



) 



2AT/Ei 



E 



'ZZ 



E' 



\zz 



R — 



where the spatial scale i?^ = lj^Ei/ma.x{t, fi^ , \g\iiBH, E^^) indicates the ex- 
citation wavelength where its anisotropy energy starts to compete with its 
exchange energy. Note that the three first constants renormalize as an ex- 
ternal field whereas E^^ constant renormalizes as an easy-axis anisotropy. 
Although the Coulomb energy Ei ~ lOOK ^ T IK in most experiments, 
the specific number: 24 = 3 x 8, which is related to the order of anisotropy 
and to the eight degree of freedom for thermal fluctuations in the case of 
SU (4) symmetry, makes the renormalization of the constant E^^ noticeable. 

As we have seen in the previous section the order parameter can be pa- 
rameterized by six or eight angles in the case z/ = 1, 3 or z/ = 2. Actually, 
not every of those rotations corresponds to a physically distinct eigen state. 
We restrict the calculation of the total bilayer energy up to a first order 
in powers of the anisotropy Hamiltonian, which means that we shall need 
only its diagonal matrix elements. But, these are real matrix elements of 
course, despite the fact that in an external magnetic field there is no time 
reversal symmetry. Hence, if the Hamiltonian is real one so real has to be 
its ground state. One generates all real eigen states from a reference state 
by rotations from the 5*0 (4) sub-group of the S'f/(4) group. This group 
has 6 parameters with two of them falling into the denominator sub-group. 
Thus, the ground state differs from the reference state by just four rotations. 
One can view locally the 8D manifold of order parameter as a composition 
of four unit vectors: magnetization of the first and the second layers and 
the two hopping-tau vectors which represent the distribution of spin-up and 
spin-down electron density over the two layers. Now the first two term in 
the Hamiltonian (|T2p are external fields acting on these four vectors. On 
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the other hand the Coulomb energy couples pairs of tau vectors via an ex- 
change interaction. This instructive picture allows us to identify only three 
special global rotations that do change the total bilayer energy. We start 
with the case v = 2 and we use a set of trial many electron wave functions 
parameterized by the three angles of relevant in our case rotations: 9± and 

n U{{}, -m{0+, 0^) cX,^c% \empty), (49) 
p 

where spins in the layer 1,2 are first rotated by 
f/|-(^,-^) = exp(z-a^)„^+ exp(-^-a^)„^, 

\ / TlT2 V / riT2 

(50) 

and then wave functions of electrons with spin ± spill over the two layers, 
the process described by two distinct angles: 9±,: 



a 



° ' a'\ ,9. /(t°-(t"\ , .9. 



i?^-(0+,^_)= J exp(z^r^).,., + J exp(-^-r^).,.,. 

(51) 

This set includes the singlet- liquid state at 9^. = 7r/2 and = and the 
canted antiferromagnetic state at 9± = 0. The order parameter reads: 

Q = URNR+U+, (52) 

with being the electron density calculated with the reference state of the 
previous section (see (0)). Now we substitute (|5^) into the anisotropic 
Hamiltonian (^7]) to find the the total anisotropy bilayer energy as: 

E^^^' = -E''cos9+cos9^ - (t + E°^) cos^9(sin 0+ + sin 0_)- 

-(/i" + E°")(cos0+ - cos^_) - l^l/is^T sin^9(cos^+ + cos0_), ^ ^ 



The minimum of this energy corresponds to three phases: a) ferromagnetic 
■i? = 7r/2, ^+ = 6*. = 0; b) spin singlet = 0, 6*+ = tt — 6'_ = 6^; and c) canted 
antiferromagnetic state otherwise, as it is shown on the Fig.3. It is identical 
to that found in the Ref . |]1T| . A line of continuous phase transitions between 
the ferromagnetic phase and the canted antiferromagnetic phase is given by 
the following equation: 



\g\y.Bli = (t + ^"^y {E^^ + H^bH) 

(54) 
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H A 




Figure 3: Phase Diagram in the Gate-symmetric case 

In the spin singlet phase the interlayer mixing phase: 9, is determined by the 
equation: 

(E'' sin ^ + t + cos ^ = (/i" + E^') sin 9. (55) 

A phase transition hne that separate the spin singlet phase from the canted 
antiferromagnetic phase is given parametrically by the equation: 

(it + sin e-E'' + (/i^ + cos e) (t + = {\g\fiBHf sin 6, 

(56) 

where 6 is determined from (0). This phase transition is a continuous one 
also. 

In the case u = 1 there is no Coulomb interaction energy and the total 
bilayer energy reads: 

^anis ^ _^gin^_^2pos^_ \g\j^^H COS ^, (57) 

The minimum of this energy is given by electron spin being directed along the 
magnetic field: = 0, whereas 6 = tan~^t/yU^. There is no phase transition 
in the case u = 1 and the only phase can be characterized as ferromagnetic 
in both the spin and the layer spaces. 
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The case u = 3 formally reduces to the case i/ = 1 although here the 
Coulomb interaction energy does not vanish identically. We find renormal- 
izations to the one-particle electron Hamiltonian whereas the total energy 
being similar to the case v = 1: 

^ams ^ _(t + 2E°^)sin^ - ifi' + 2E^')cose - \g\i2BH cos^, (58) 

There is no phase transition in this case either. 



Anisotropic energy gap of one Skyrmion 



In this section we find an anisotropic part of the total skyrmion gap energy. 
The non-homogeneous order parameter that represents one skyrmion is given 
by the Belavin-Polyakov (BP) skyrmion solution for |Q| = 1 0]: 



Qbp{zz) 
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zR 
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2 


[ zR 


R' 



(59) 



where only the shown above four matrix elements differs from those in the 



electron density matrix (31). The skyrmion order parameter has to be 
rotated by a homogeneous matrix RU calculated in the previous section in 
such a way that the order parameter far away from the skyrmion center 
minimizes the anisotropy energy. In addition to this rotation we have to 
allow all homogeneous rotations from the denominator sub-group S that 
actually transform the BP skyrmion order parameter (|59|): W. Thus, a 
general skyrmion order parameter reads: 

Q(f) = RUWQbp{zz)W+U+ R+ . (60) 

2. Here the matrix W is parameterized by 



First, we consider the case v 
seven angles: 



/ cos/?/e*(^^+°^) 
-sin/?/e^(-^-f+"^) 





sin;g/e*(^/-"-f) 
cos/3/e*(-^/-"/) 










COS /5ee*('''=+"^) 
- sin /?ee^(->+°'=) 



\ 



sin/^ee*^^-"'^) 
cos/?ee*(-^^-°^^ / 
(61) 
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The additional seventh parameter angle of the denominator group just rotates 
the coordinates: z — > e*''"'^. We find by explicit calculation that the skyrmion 
anisotropic energy does not depend on parameters 7e, 7/ and 77 whereas 
= and = vr correspond to the skyrmion energy minimum. Thus we 
shall express the anisotropic skyrmion energy in terms of the two relevant 
parameters /9e and /9j that rotates the core of skyrmion in the empty-electron 
space and fiUed-electron space correspondingly. Besides these angles the 
skyrmion energy depends on free parameters entering the BP solution. The 
calculation here has been performed only for skyrmion topological index: 
|Q| = 1 ( p9D and in this case there is only one such parameter, namely, the 
radius R of the core of skyrmion. There are few different spatial integrals 
that we encounter in calculation. And only one of them is logarithmically 
divergent. We calculate the skyrmion anisotropy energy with the logarithmic 
accuracy. It means that parts of the anisotropy energy coming from different 
anisotropy sources are all multiplied by the same spatial integral, which we 
denote by a constant A: 

/ p \ 2 p* 
A-{-yo,-. (62) 



where B* is the inverse mass of Goldstone excitations in the model (|35|j53| ). 
Also, the Zeeman energy of skyrmion reads: 

-^skyr _ _ j\\^g\^^^]^ ( sin'(9(cos/3/ — COS /3e) (cos 6'+ — cos6'_) + 

9 -\- 9 \ 

+ sin'i9(cos 6'+ + cos 9_) — cos'i9(sin/5/ + sin /3e) sin ^ ^ j 

(63) 

The hopping energy of skyrmion reads: 

^skyr ^ ^ ^ ^Ox) }_ cos^9(cOS/?/ - COS /?e) (siu ^+ - sin^_) + 

9 9 \ 

+ cos -(9 (sin 6*+ + sin ) — sin -i? (sin /5/ + sin /?e) cos ^ j 

(64) 

The gate asymmetry energy of skyrmion reads: 

E'^y' = A8g = A{^i, + E^') {{cos 9+ - cos0_)- 

-(cos/?/ — cos/3e)(cos6'+ + cos6'_ ^ ^ 

2 
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And finally the anisotropic Coulomb energy of skyrmion reads: 

^skyr ^ ^g^^ ^ ^^zz 1(1 + cOS/3/COS/3e)(l + COS 0+ COS0_) + 

^ 2 ^ X (66) 

+2 cos cos 9_ — - sin /?j sin (3^ sin 0+ sin 9_ \ . 

There is also a contribution to the skyrmion energy coming from an non- 
homogeneous BP electric charge distribution inside the skyrmion core 



1 /■ /■ R2 ^2 E>2 o„2 „2 

^ 27 i 7r(r2 + i?2)2^|^-._ -,|^(^/2 + ^2)2 64 ^ ' 

The minimum of the total anisotropic skyrmion energy (|63| . |6^ , |65| , |66|) : S^^^^ = 
£z + ^h + £g + £zz, over the two parameters /3e and /?/ was found numerically 
and is denoted as: £^in ■ Next, we find a minimum the total skyrmion energy 
including (0): E'^^y'^ = AS^-^^ + E^'^\ with respect to the skyrmion radius 



This formula holds in the limit S''^'^'^^' <^ Ei and, thus, the second term is much 
smaller then the first term in the Eq.(p8|) as they are calculated. But, it is 
important in the case of antiskyrmion Q = — the gap is only a relatively 
small anisotropic energy. The resulting anisotropic part of a skyrmion gap is 
shown on the Fig. 4 in the case /x^ + E^'^ = 0. Note the two prominent cusp- 
like lines in the skyrmion gap sheet coincide with the two phase transition 
lines from the Fig.3. A view of another cross-section of the skyrmion gap: 
\g\fiBH = 0, is shown on the Fig. 5. A minima here also coincide with the 
canted-antiferromagnetic phase. A skyrmion in the ferromagnetic state is a 
spin-skyrmion with the spin rotations being localized inside one of the two 
layer, whereas a skyrmion in the spin-singlet state is a layer-skyrmion with 
the electron density being distributed over the two layers. 

In the experimental setup ^ gskyr formulas can be com- 

pared with the experimental results only qualitatively. But they found a 
profound disappearence of the thermal activation gap maximum at some in- 
terval on the z/ = 2 line. In our theory this would correspond to a minimum 
in the skyrmion gap and we suggest that this indicate the canted antiferro- 
magnetic phase. 
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Gap 




Figure 5: Anisotropic part of Skyrmion Energy in the case of vanishing 
Zeeman energy. Spin-singlet phase is on the right side and the canted anti- 
ferromagnetic phase is on the left side 
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In the case z/ = 1 we parameterize general rotations from the denominator 
sub-group by four angles in such a way that the electron density matrix 
becomes: 



WQbp{zz)W' 



1 f 


z 


2 


Rz\f) 


B? + 


z 


2 V ^^(/l 


RV){f\ 



where 



I/) 



cos — , sm — cos ae , sm — sm ae 
2' 2 '2 



(69) 



(70) 



Once again we find that the skyrmion energy does not depend on the param- 
eters Ai and A2. A straightforward calculation shows then that the Zeeman 
skyrmion energy is 



A = 2A\g\^BH y- ~ sin^ ^ cos^ a 



(71) 



the Hopping skyrmion energy is 

At — sin (3 sin a + sin 6 



^skyr 



1 2 • 2 

1 + COS a sm — 



whereas the gate asymmetry skyrmion energy is: 



E. 



skyr 
G 



— cos l3 



1 • 2 2 

1 — sm — cos a 



e 



cos a sin - + cos 9 



(72) 



(73) 



Note that there is no Coulomb energy in the case v = 1. Searching for minima 
of Eg. (|71| , [72|J73D varying parameters j3 and a and fixing 9 = tan~^t//i^ we 
find the minimal anisotropy energy of skyrmion to be: 

A = ^^^E^ + min (2^^^ + ^2 ^ 2\g\fiBH'^ • (74) 



The skyrmion gap is given by the same formula as in the case u = 2: 
Generally, there is no prominent minima in the skyrmion gap here. Nevether- 
less in the limit t <^ \g\fiBH such a minimum occurs. 

The case u = 3 reduces to the case u = 1 ii one simply to renormalize the 
constant in (^ij) as it was explained in the previous section. 
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Conductivity Activation Energy 



The conductivity activation energy was measured in experiments for bilayer 



system p] and in one layer system [|14[. It was known for quite some time 



that this energy is considerably less than a typical exchange constant e^/ kIh- 
This fact was the main motivation for the experimental search of topological 
excitations. If one considers an act of creation of skyrmion and anti-skyrmion 
pair at large separation then one easily gets the pair energy: Ei = 1/2Eq. 
This is definitely less than the creation of electron-hole pair at large separa- 
tion: Eq, but still of the order of / kIh-, for one layer. In a bilayer system 
the prominent minimum in the minimal activation energy for skyrmion anti- 
skyrmion pair is still of the order of in spite of the fact that experimental 
conductivity activation energy goes to zero very sharply 

This controversy can be overcome by considering the creation of neutral 
antyskyrmion and electron pair at large separation. First, we consider a one 
layer case. The energy of the additional electron with reversed spin does 
not contain exchange Coulomb energy and consist from anisotropy energy 
only. The energy of anti-skyrmion also has only anisotropic energy (p^). 
The anisotropic energy of electron can be neglected in the limit of large 
ratio: Ei/\g\^BH , which holds in most experiments. The total anti-skyrmion 
energy is: 

72 -* -I 2 

\g\fiBB /(l-cos/?(r))^^ + - / — ^— curl n'(r] curl n'(r')(f rd^f*. 
J 2'Kljj 2 J K\r — r I 

(75) 

. Plugging the BP solution we calculate this energy to be: 

2|,|,,B«Mo4 . (76) 

Minimizing it further with respect to radius of the anti-skyrmion core R we 
find the core radius: 

R=h ( V' (77) 

2''\32KlH\g\^^Blog{R*/lH)J ^ ' 



and the activation energy: 



1/3 



^ = '^-\\g\^JiBH\og{W/lu)[^^] 1 . (78) 
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The upper limit under the logarithm R* is defined by the validity of BP 
solution. Essentially, it is a distance where gradient energy becomes of the 
order of Zeeman energy: 



Therefore the logarithmic factor is of the order of unity. If it were rather large 
one would need to compare the energy of Q = — 1 and Q = —2 antiskyrmions, 
with the latter being logarithmic free. In any case an important point is that 
the activation energy is magnetic field dependent: 

A = KH^/^, (80) 



which in agreement with the various experimental results |TJ] . The absolute 
value of the constant K (|7Spin this relation also conforms experiment [|T^ . 
Note that the Zeeman energy of electron is much less that A (|80D . In the 
case of double layer system the situation is more complicated due to large 
number of parameters defining the if^™*. Nevertheless, assuming anisotropy 
energy is small and neglecting for the same reason the anisotropy energy of 
an electron we get the deep minimum in the canted antiferromagnetic phase 
in accordance with the experimental result @]. It should be noted that in the 
process of creation of anti-skyrmion electron pair the total topological charge 
of the system is not conserved as opposed to the skyrmion antiskyrmion pair. 
Therefore, this process goes tentatively on the sample boundary. Also the ex- 
istence of magnetic field makes it possible to violate electron-hole symmetry 
usually assumed in theories used the basis of projected on the lowest Landau 
Level wave functions. The violation of the electron-hole symmetry is related 
to the topological charge and gives rise to the skyrmion vs antiskyrmion 
energy difference due to the topological term in the action (p5|). 
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